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Abstract: As a first step towards the application of the antenna subtraction formalism to 
NNLO QCD reactions with massive quarks, we determine the real radiation antenna func- 
tion and its integrated counterpart for reactions of the type S — )• QQqq, where S denotes 
an uncolored initial state and Q, q a massive and massless quark, respectively. We compute 
the corresponding integrated antenna function in terms of harmonic polylogarithms. As an 
application and check of our results we calculate the contribution proportional to a'^eqNf 
to the inclusive heavy-quark pair production cross section in e~^e~ annihilation. 
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1 Introduction 

The calculation of differential cross sections and distributions in perturbative QCD beyond 
the leading order requires, apart from the renormalization of the QCD coupling and the 
quark masses, methods to regularize and handle the infrared (IR) divergencies that appear 
in the intermediate steps of such computations. One general approach is to construct 
subtraction terms such that, after adding/subtracting these terms, the IR singularities 
are regulated and cancelled in tree-amplitudes involving the radiation of real massless 
partons and in associated loop amplitudes in the computation of IR safe observables (up 
to factorization of collinear initial-state singularities). 

For calculations at NLO QCD a widely used version of this approach is the dipole 
subtraction method for massless QCD [1] and for QCD with massive quarks and other 
colored massive particles [2-4], which was slightly modified in [5-9] and has found a number 
of computer implementations [7, 9-12]. Other NLO subtraction methods were constructed 
and applied, too, including those of [13-16], and the antenna method (see below). 

Computations of differential cross sections at NNLO QCD involve three types of con- 
tributions: squares of tree-level double real emission amplitudes (with n + 2 final-state 
partons), interferences of one-loop and tree-level amplitudes {n + 1 final-state partons) and 
of Born, one-loop, and two-loop amplitudes (n final-state partons). For general discus- 
sions of the IR structure, see [17-19]. Various techniques have been devised to handle the 
IR divergences of these individual contributions. These include the sector decomposition 
algorithm [20-22], the antenna formalism [23-25], and the subtraction methods [26-34]. 
Application to reactions at NNLO QCD include pp ^ H + X [32, pp ^ W + X 
[36, 37], e+e" — >• 2 jets [38, 39], and e"'"e~ — )■ 3 jets [40-43], where the jet calculations just 
mentioned were made for massless partons. 
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The antenna method [23-25], which we employ in this paper, was first worked out fully 
to NNLO QCD for e~^e~ annihilation into massless final-state partons. The general set-up 
of this approach applies also to colored initial states and/or massive colored particles in the 
final state. The extension to processes with initial-state and massless final-state partons 
at NLO QCD was made in [23, 24, 44]. For reactions with initial-state and massless final- 
state partons at NNLO QCD, results were presented in [45-47]. For processes with massive 
quarks Q in the final state, the antenna subtraction terms at NLO QCD were explicitly 
worked out for colorless initial states and for the hadronic reactions /11/12 — )• QQ, QQ+ jet 
in [48] and in [49], respectively. As is well-known, the IR singularity structure of the 
matrix elements for reactions with massive colored particles is less entangled than that of 
their massless counterparts but, on the other hand, the (analytical) computation of the 
integrated subtraction terms is more difficult. 

In this paper we are concerned with the production of a heavy-quark pair QQ by an 
uncolored initial state S at NNLO QCD, i.e., we consider reactions of the type 

S^QQ +X, (LI) 

at order a^. This includes the production of a pair of heavy quarks by electron-positron 
annihilation, e'^e~ — )• 7*,Z* — )• QQX, and the decay of a color and electrically neutral 
massive boson of any spin into QQX. The following ingredients are neccessary for the 
computation of arbitrary differential distributions to order a^: i) The amplitudes S — )• QQ 
to order a^. They are known in analytical form for S = vector [50, 51], axial vector 
[52, 53], scalar and pseudoscalar [54]. ii) The tree- and one-loop amplitudes for S — )■ QQg. 
The one-loop amplitudes can be computed with standard methods and are known for 
e^e~ annihilation, i.e. S = j* , Z* [55-57]. iii) The tree-level amplitudes S — QQQQ, 
QQgg, and QQqq, where q denotes a massless quark. Apart from the tree-level amplitudes 
S —7- QQ and S — t- QQQQ, the matrix elements give rise to IR singularities, which are 
regulated within the above-mentioned subtraction methods by appropriate counterterms. 
The calculation of the differential cross section dcTNLO for S decaying into two massive 
quark jets is standard. The contribution of order to the two-jet cross section is given 
schematically, using the notation of [39], by 



Here c^^^j^lo ' ^"^nnlo ' ^"^nnlo denote the contributions from the tree- level amplitudes 
QQgg and QQqq (and QQQQ, which does not require a subtraction term), the one-loop 
QQg amplitude, and the two- loop QQ amplitude, respectively. The term c^ctnnlo ('^'^nnlo) 
is a subtraction term that coincides with da^j^^^Q (d^js^J^Lo ) singular limits. 

The IR singularities of the two-loop term dcNNLO explicitly known within dimen- 
sional regularization [50-54]. The construction of dcr^NLO da^^'^Q depends on the 
subtraction method used - the integration of these terms over the four- and three-parton 
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phase spaces and d<I>3, respectively, is in any case a difficult task. To our knowledge 
this has not yet been done for massive quarks in analytical form. As mentioned above, 
we shall use the antenna framework. As a first step in the computation of da{QQ) to 
NNLO QCD within this approach, we determine in this paper the subtraction term for 
the QQqq final state and its integral over the four-parton phase space. The new aspect 
of this computation is the analytic integration of the massive tree-level antenna function 
associated with the process 

l*{q) ^ Q{pi) Q{P2) + q{P3) q{pA) • (1-3) 

over the full four-particle phase space. This (integrated) antenna function is not only of 
relevance for the specific process at hand, but serves also as a building block for constructing 
subtraction terms for other processes (1.1) within the antenna formalism. 

The paper is organized as follows. In Section 2 we determine the antenna function for 
(1.3) and in Section 3 we integrate this function analytically over the four-parton phase 
space. As an application and check of our results, we compute in Section 4 the cross 
section for the inclusive production of a massive QQ pair plus Nf massless quarks by e'^e~ 
annihilation via a virtual photon - more precisely, the contribution of order CqolNf to 
this cross section. Section 5 contains a summary and outlook. 



2 Antenna subtraction at NNLO QCD 

In the following, we restrict our attention to the case of reaction (1.3) where, for definiteness, 
we consider one massless quark flavor q in the flnal state. As mentioned above, we focus 
on constructing a subtraction term which coincides with the squared matrix element of 
7* — )• QQqq in all single and double unresolved limits. In fig. 1 the Feynman diagrams are 
shown that are associated with this process at order a^. The corresponding contribution 
to the cross section for 2-jet production may be written as follows: 



da 



R,QQgq 
NNLO 



ATraiAnasY [N^ - l) d^4{pi,p2,P3,Pi;q) A'^\pi,P2,P3,P4. 



X <4 



M° - 
\QQ<i 



where the matrix elements M^^^-q and M^^^q- correspond to the diagrams Ci, C2 and 
C3, C4, respectively. The dependence on the electromagnetic and strong coupling and 
the dependence on the number of colors Nc are extracted from the matrix elements; eq 
(cq) is the electric charge of the massive (massless) quark in units of the positron charge. 
Summation over all spins is understood. In the formulae below, the polarizations of 7* are 
summed, but not averaged. 

The phase space measure in d = 4 — 2e dimensions is 



d^4{pi,P2,P3,P'i;q) = M 



12-3d 



n (27r)^-i2p' 



{271^6^'^'^ 



4 

E 



Pi 



(2.2) 
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Figure 1. Feynman diagrams contributing to 7* — ?> QQqq at tree-level. Bold (thin) lines refer to 
massive (massless) quarks. 



where ;U is a mass scale. The jet function Jm^ in (2.1) ensures that only configurations are 

taken into account where n outgoing partons form m jets. 

2 



Because the squared matrix element 



- 

qQQq 



does not involve infrared singular con- 



figurations, no subtraction is required. Its contribution to the QQ production cross section 

is given in [581. The same holds for the interference terms between A^9, and A^''^/=,-. 

° ^ J QqqQ qQQq 

Moreover, due to Furry's theorem these terms yield a vanishing contribution to the cross 
section if the observable under consideration does not distinguish between quarks and 
antiquarks. 

In the framework of antenna subtraction the main building blocks for constructing 
NNLO subtraction terms are the antenna functions, which can be derived from physical 
color-ordered squared matrix elements for tree-level 1—7-3 and 1—7-4 processes and one- 
loop 1 — >• 3 processes. A detailed and completely general analysis of how the subtraction 
terms are constructed from the various antenna functions is given in [25] for massless final 
state partons. This procedure applies also to the case of massive quarks. Its application 
to the specific process 7* — >• QQqq is outlined below. The three-parton tree-level antenna 
functions with massive quarks have been calculated in [48, 49], whereas the four-parton 
tree-level and three-parton one-loop antenna functions that involve massive quarks are still 
missing. 

In the case of 7* — >• QQqq the construction of the subtraction terms may be divided 
into two parts. In a first step, the single unresolved configurations are subtracted. Within 
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the antenna method, the corresponding subtraction term reads 



da 



S,a 

NNLO 



Aira (Airasf Cq [N^ - l) d$4(pi,P2,P3,P4; 

ii;3^lQ, 3,, 4,-) AO ((T3)q, (43)^, 2q) (^^3,]5r3,P2) 
+ S3°(2q3„45) ^"(Iq, (34)^, (24)q) jf 



(2.3) 



The massive quark-antiquark antenna function ^{iq, kg,jQj and the quark-gluon antenna 
E^{iQ,jg, kg) with a massive radiator quark are given in [48]. The momenta and pj^ are 
redefined on-sheh momenta, constructed from hnear combinations of the momenta pi, pj 
and pk- The tree-level two-parton matrix element squared (summed over colors and spins, 
with the photon coupling and Nc factored out) is 



A^qq(7* ^QQ) =4 [(1 - 6) g2 + 2m^] , 



(2.4) 



where m denotes the mass of Q. In a second step, the double unresolved configuration, 
where both q and q become soft, has to be subtracted. In the case at hand, the appropriately 

2 



normalized squared matrix element 



- 

QqqQ 



can be used as subtraction term. In the 



terminology of [25] this is the antenna function B2{1q, 3^,4^, 2q) associated with the color- 
ordering QqqQ, where a color-connected massless quark-antiquark pair is radiated between 
a pair of massive quarks. More precisely, this antenna function is defined by 

2 



^4(lQ)3g,4g,2Q) 



- 

QqqQ 



QQ 



2 ' 



(2.5) 



where the normalization factor is given in (2.4). Obviously (2.5) has the appropriate 
behaviour in the singular double unresolved limit where q and q become simultaneously 
soft. However (2.5) contains also singularities due to single unresolved limits, which have 
to be subtracted from the antenna function. In the end, the corresponding subtraction 
term for the double unresolved configuration reads 

2 



'^^NNLO = 4™ {^T^asY eg {N^ - l) d^i{pi,P2,Pz,PA; Q) 



-^3(2q3„4,-) A[|(1q,(34)^,(M)q 



where piki and pjki are linear combinations of the momenta pi, pj, pk and pi. 
The subtracted differential cross section 



(2.6) 



, R,QQqq _ , S,a 
""NNLO ""NNLO 



da 



S,b 

NNLO 



(2.7) 



is free of IR divergences and can be integrated over the four-parton phase space numerically 
in d = 4 dimensions. 
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For the antenna function we find 

-B4(lQ,3g,%,2Q) = + 2m2) { ^^^^ [*12S13 + S12S14 + S13S23 + S14S24] 

H [S12S23 + S12S24 + S13S23 + S14S24] 

■534'S234 

1 r 

+ 2 2 |.2si2Si3Si4 + S13S14S24 + S13S14S23 
•534'5l34 

- S?3'S24 - S?4'S23] 

1 r 

+ 2 2 I.2S12S23S24 + Sl3S23'S24 + S14S23S24 
•534'5234 

- 513524 - Susls] 

H [2s?2 + S12S23 + S12S24 + S12S13 + S12S14] 

S34S134S234 

+ -o [ - •Sl3si4 - S14S23 - S13S24 - S14S23 

•S34'Sl34'S234 

+ S13S14S23 + S13S14S24 + S13S23S24 + S14S23S24 

2S12 



2S12S13S24 - 2si2Si4S23 + 



•Sl34'S234 



2 



8S13S14 , 8S23S24 



+ m — K 1- 



q2 --i2 q2 q'2i 

*34*134 '^34*234 *134 *234 



2 2 

[Sl2 + S23 + •S24] n— [Sl2 + Sl3 + S14] 



S34S134 ■534S234 

2 

H [4si2 - Sl3 - Sl4 - S24 " S23] 

■S34'Sl34'S234 

8 

2 [S14S23 + S13S24J 

•S34'Sl34'S234 

- + ] + ^(^) ' (2-8) 

Vs34Sf34 S34S234/ J 

where Sij = 2pi ■ pj and Sijk = Sij + Sik + sjk- For the sake of brevity we have not 
written down in (2.8) the terms of order e. For the numerical computation of (2.7) in 
d = 4 dimensions, only the four-dimensional antenna function is required. However, 
the integrated antenna function, which we compute in the next section, must be determined 
with i?4 in d dimensions. 

For completeness, we derive also the behaviour of the antenna function B^ in the 
double soft limit ps — )• Xps, p^ — )• Ap4 with A — )• 0. We obtain B^ = A~^S'i2 + 0{X^^) with 

2 

Sui^gAg) = -Q—. TT ; T (S12S34 - S13S24 - S14S23) 

S§4 (Sl3 + S14) (S23 + S24) 

2 /S13S14 - S34m2 S23S24 - 534™-^ ' 



■534 V (Sl3 + Suf (S23 + S24:f J 

The double-soft factor (2.9) is not identical to the respective factor S^"^^ for a massless 
quark Q given in [25], but agreees with this factor in the limit m — t- 0. A respective 
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difference between the m ^ versus m = case was sliown in [34] for QQgg final states in 
the double-soft gluon limit. 



3 The integrated QqqQ antenna function 

The introduction of subtraction terms must be counterbalanced by adding their integrated 
counterparts, cf. (1-2). In the context of the antenna method this implies the analytic 
integration of the antenna function i?4(lQ, 3g, 4g, 2q) over the four-parton antenna phase 
space d^Xg^i^Q associated with a massive and a massless quark-antiquark pair. The corre- 
sponding integrated antenna function ^Q^qQ is defined by 

^Q.,q(9''"^'^'^) = (8^' (4vr)-^e^^-)' J d<Px^^^^ B',{1q,3„A„2q) . (3.1) 

The antenna phase space d^Xg -q is closely related to the full four-particle phase space 
(2.2): 

d^i{pi,P2,P3,Pi,q) = P2{q^,rn)d^xQ^^Q, (3.2) 

where P2{q'^,m) denotes the integrated phase space for two particles of equal mass m in 
d = 4 — 2e space-time dimensions: 

The integration of the antenna function over the antenna phase space has to be performed 
analytically in d = 4 — 2e dimensions. 

3.1 Reduction to master integrals 

In order to calculate the integrated antenna function (3.1), we first reduce the number of 
integrals to be computed by exploiting linear dependences between phase-space integrals 
with the help of integration- by-parts (IBP) identities [59]. These identities were originally 
derived for loop integrals, but are also very useful in computing phase space integrals 
[60, 61]. 

Using the Cutkosky rules [62] we introduce two massive and two massless cut-propagators 

^ = 27ri5+{p^-m^) = - 1^^^ - 1^-^, i = 1, 2, (3.4) 

Di Pi — ^ + *0 Pi — — lO 

= 27ri5+ (p2) = - i = 3, 4. (3.5) 

Then we can write: 

#4 {Pl, P2, P3, P4, q) = ■4^27T)^d '^^"'^ " Z^^^^^ 11 ^^'^^ 

Next we express each of the invariants Sij, Skim in the d-dimensional version of (2.8) in 
terms of the four cut-propagators and five further appropriately chosen propagators and 
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scalar products, such that each term of B^^-q can be viewed as a four-particle cut through 
a three-loop vacuum polarization diagram, which may contain irreducible scalar products 
in the numerator. It was shown in [60] that IBP reduction of such integrals can be carried 
out in the same way as for loop integrals, using that integrals where at least one of the four 
cut-propagators is missing in the integrand do not contribute to the original cut-integral. 
Using the implementation AIR [63] of the Laporta reduction algorithm [64] we decompose 
the integral (3.1) along these lines. As a result we can express the integrated antenna 
function in terms of five independent integrals (master integrals): 



(87r2 {AT:)-'e'^^Y 



64 + 208z + 48^2 - lOz^ + 1 



z{l-zf 



+ 



2 (32 - 448z - 112^2 + 133^^ 



49z^ 



+- 



3z(l -z)3 
2 ( 1720 + 5656z - 164^^ 



1 

e 

163z^ + 2142^ 



+ 



9z(l - zf 
4 (48 + 56z + 6^2 - 5^3) l 



+ 



+ 



z{l - zf e2 

4 (72 + 56z + 234^2 - lOlz^) 1 
3z(l - z)3 7 

8 (1296 + 1222z + 153^2 - 241^^) 
9z(l - zf 

4 (16 + 64z -I- 26^2 - z'^) 1 
z{l-zf ^ 

4 (8 - 232z - 50^2 + 13^3) i 
2,z{l-zf e 

8 (440 + 1730z + 247^2 + 132^^ 



+ Oie) 



9z(l - z) 
{qY"T,{q\m\e) 



+ 0{e) 



+ 



8(4-^2)1 4 (20 + 12z -h 49^2) 



+ 0{e) 



3(1 - z)z e 

2 ™2 



- + 



9(1 - z)z 



+ Oie) 



X Ti{q^,m^,e 
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+ 



2(4-22)1 4(8 + 92 + ^2) 



3(1 -z) e 



9(1 -z) 



+ 0{e) 



(3.7) 



where we introduced the dimensionless variable z = Am /q . The five master integrals are 



Ti{q^,m^,e) = ^^^^^ ^ I d^4{Pi,P2,P3,P4;t 



(3.8) 



T2{q ,m ,e 



T^{q^,rr?,e) 




Sim 



d^i{puP2,Pz,Pi]q) si3 



d^i{pi,P2,Pz,Pi]q) S134 



(3.9) 



(3.10) 



Ti{q^,w?,e 




d^i{pi,P2,P3,P4.]i 



Sl34 



(3.11) 



T5(g^m^e) 




d<^4ipi,P2,P3,PA;< 



•S134S234 



(3.12) 



In these diagrammatic representations bold (thin) lines refer to massive (massless) prop- 
agators. In the case of T2 and T3, the invariants to the left of the cut-diagrams denote 
numerator factors. Note that Ti is just the d-dimensional phase space volume associated 
with two massless and two massive (equal mass) particles. The five integrals are all finite. 

We used also the package FIRE [65] for an independent check of the above reduction. 
Equation (3.7) shows that the integrals Ti, T2, and T3 have to be computed to order e^, 
while it is sufficient to compute T4 and T5 to order e. 

3.2 Analytic computation of the master integrals 

In the integrands of Ti, T2 and T3 there are no denominators present, and the respec- 
tive numerator factors depend just on a subset of phase space momenta. Based on this 
observation, we first rewrite the phase space (i$4 in terms of the convolution formula 



1 

d^A{pi,P2,P3,P4;q) = TT dM"^ d<^2{P4,k;q)d^3{pi,p2,P3;k), (3.13) 
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where k"^ = M^. Using this relation along with standard identities and integral represen- 
tations of hyper geometric functions [66], we find that the first three master integrals can 
be expressed in terms of hyper geometric functions 3F2: 

Ti-{q)\^-^) |2 vr r(3 - 3e) r(4 - 4^) 

X 3F2 (-^ + €,-2 + 3e, -3 + 4e; e, -1 + 2e; z 



)-12+8e ^-5+3e 



r(l-e)3r(-l + e) 
r(3 - 3e) r(2 - 2e) 



X zi-Si^2( ^,-l + 2e,-2 + 3e;2-e,e;z 



,-15+lOe -5+3e - e) r(-l + e) 



r(2 - 2e) 

xz'-'^3i^2(^-e,e,-l + 2e;3-2e,2-e;z) }, (3.14) 



T2 = {q^f /"/^y j2-12+6e^-5+3e r(l e)^ 



IT 



) \ 3r(3 -3e)r(4-4e) 

X 3F2 f-^ + e, -3 + 3e, -4 + 4e; e, -2 + 2e; z 



-13+8e _-5+3e 



r(l -e)3r(-l + e) 
3r(3 - 3e)r(2 - 2e) 



X z^-^ 3F2 ( ^, -2 + 2e, -3 + 3e; 2 - e, -1 + e; z 

, ^-16+ioe ^-5+3e r(2 - g) r(-2 + e) r(-i + g) 

+ ' " f(232^^ 

x/-2Si^2(^-e,e,-l + 2e;4-2e,3-e;z) }, (3.15) 

'~^'^W) I r(3-36)r(4-4e) 
X 3i^2 (^-^ + €,-2 + 3e, -4 + 4e; e, -1 + 2e; 2; 

^ 2-11+8. ^-5+3. r(i-6)3r(-i + 6) 



3r(3 - 3e)r(2 - 2e) 
X zi-^gFaf ^,-l + 2e,-3 + 3e;2-e,e;z 
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+ 



2-15+10e^-5+3. 



. r(i-6)r(-i + 6)^ 

r(2 - 2e) 



X z 



.2-2e 




-e,e,-2 + 2e;3-2e,2-e;z 




(3.16) 



We use the package HypExp2 [67] for expanding the hyper geometric functions in e to the 
required orders. For the above three master integrals the result of these expansions is given 
in Appendix A in terms of harmonic polylogarithms (HPL) [68, 69]. In [34] the phase space 
volume T\ was also computed by expansion in e in terms of harmonic polylogarithms. 

The remaining two master integrals T4, T5 are computed using the differential equa- 
tions method [70-72]. (For a review see [73].) By differentiating T4, T5 with respect to 
z = Am? /q^ we obtain linear combinations of integrals which in turn can be reduced by 
IBP identities to the above five master integrals under consideration. In this way we derive 
inhomogeneous first order differential equations in z for T4 and T5. 

The inhomogeneous part of the differential equation for T4 only depends on the three 
master integrals Ti , T2 , T3 whose series expansions in e were already derived above in terms 
of HPL. We expand this equation for T4 in e: 



In this way we obtain first order differential equations for the coefficients and , 
whose inhomogeneous parts are determined in terms of HPL. In each case the general 
solution is composed of the general solution of the homogeneous equation, which contains 
a constant of integration, and an integral over the inhomogeneous part. The differential 
equations for are obtained and solved in a similar way. Here the inhomogeneous part 
depends on the other four master integrals. 

The integration over the inhomogeneous parts is carried out using the package HPL [74] . 
After partial fraction decomposition and expansion of shuffie products we can write the 
respective integrand such that all terms containing HPL are of the form kf{zy H{...,z), 
where A: is a constant, f{z) is one of the functions j^, \-, appearing to a positive 
integer power j, and where H{...,z) is a HPL of weight w and argument z. In the case 
of j = 1 the primitive is k times an HPL of weight w + 1, which follows directly from the 
definition of the HPL [68]. In the remaining cases where j 7^ 1 we can perform partial 
integration and partial fractioning in sequence until each term is either of the above form 
with J = 1 or is just an algebraic function of z. 

In order to fix the constants of integration of the integral T5, one can consider the 
massless limit z = where T5 remains finite. In this limit the integral was computed in 
[61]. We use that result as a boundary condition in order to determine the constants of 
integration for T5. 

In the case of the differential equations obtained for T^^^ and T^^\ the massless limit 
z = and the threshold limit z = 1 can not be used as boundary conditions for the 
determination of the integration constants. Instead we choose an appropriate integral 
which is known to vanish in the limit z = 1 and which can be expressed in terms of the five 




(3.17) 
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master integrals by IBP reduction. Via this reduction we use the latter limit as boundary 
condition for fixing these integration constants. 

For T4 and the result of their expansion to oder e in terms of harmonic polyloga- 
rithms is also given in Appendix A. 

For the integrals T4 and T5 we performed numerical cross checks using VEGAS [75]. A 
strong analytical check of all five master integrals is provided by analysis described in the 
subsequent section. 

In the following we use the variable 



1 - 

y^—^^j=. (3.18) 

Inserting our results for the master integrals into (3.7) we obtain as our main result the 
integrated antenna function expanded in e: 



/^2x 2. 



-i I - 23 13 l + 2y 

+ 7 1 V 9 ^ 36(1 -y)~ 36(1 + y) ^ 6 (1 + 4y + 2/2) ) 

^^(0,0; 2/) 



3 3(1-2/) 3(1 + 2/) 
2 2 2 



3 3(1-2/) 3(1 + 2/) 



43 27r2 27r2 27r2 y 

+ 7T7^ r + 



36 9 9(1 - 2/) 9(1 + 2/) I + Ay + y'^ 
883 - 157r2 757 - 157r2 271 - 157r2 2 + 72/ 



108 108(1 - 2/) 108(1 + 2/) 3 (1 + 42/ + 2/2)' 



18 (1 + 42/ + 2/2) 7 V9 (l-2/)2 36(1-2/) 

35 5 + 222/ 



36(1 + 2/) 6 (1 + 42/ + 2/2) 
57r2 57r2 57r2 



F(0,0;2/) 

^(-i;y) 



3 3(1-2/) 3(1 + 2/) 

86 - 57r2 57r2 57r2 82/ , 

9 + 9(1^ + 9(rT^+ 1 + 42/ + 2/2,' ^ '^^ 
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- + 26 _ 4(1 + 2y) \ 

,9 9(l-y)^9(l + y) 3il + 4y + y^)) 

32 

i/(3;y)--if(l,l;y) 



3 3(1-2/) 3(l + y)y ^ '^^ 3 

20 20 20 . ^ ^ , 

H{-2,0;y) 



3 3(1 -y) 3(1 + y) 
9 9(1 -y) 9(1 + 2/); ^ '^^ 

/'32 32 32 \ rr^ . ^ ^ 

58 13 47 1 + 22/ \ 

y + 18(1 - 2/) " 18(1 + 2/) " 3 (1 + 42/ + 2/2) J ^ ' '^^ 

16 16 16 \ „ , 

, 22 22 22 \ ^ , 

/ "2 32 32 \ ^ , 

/28 28 28 \ , , ^ , 

+ 't- 3(131;)- 3(171;)) '^(-i' 0^") 
-0-t4i;-4i;)«(-1'1-°^^) 



0-T^-4i;)«(i.-i.''^«) 



2 2 2 



3 3(1-2/) 3(1 + 2/) 



^(1,0,0; 2/) 



707 - 1137r2 + 468 C(3) ^ 797r2 + 468C(3) 777r2 - 468 C(3) 



108 108(1 +2/) 108(1 - 2/) 



2(1 + 42/) 12 - TT^ - 362/ ~ 27r^y ^ 
(1 + 42/ + 2/2)2 6 (1 + 42/ + 2/2) 



(3.19) 
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4 The correction of OL^e^Nf to the ratio R 

As an application and check of our results of Section 3, we consider the ratio 

to order and to lowest order in a. Here cjpt = e^/(127rg^) is the massless Born cross 
section for e^e~ — )• 7* — )• In the following, we consider one heavy quark, carrying 

the electric charge eg, and Nf massless quark flavors. Here we are only interested in 
the contribution proportional to a'^CgNj to the ratio (4.1). This contribution is gauge- 
invariant and IR finite. Apart from the tree-level contributions, which are closely related 
to the integrated antenna function of Section 3, this term receives a two- loop contribution 
which was computed within dimensional regularization first in [50]. Using this result and 
our result of Section 3, we can check the IR poles of B^^-q. Furthermore we compare 
this contribution to R with the previous result of [76], which was obtained in d = 4 using 
different methods. 

Throughout this section we use the subscripts a'^CgNj or a'^Nf when referring to the 
contribution of these terms to a given quantity. We have 



aie^e- ^ ^ QQ + X),.,. = ^ ^ ^'^^ E , (4.2) 

with the lepton tensor 

Lf^u = 4 {ki^k2y + kiyk2^_i - g^,u ki ■ k2) , (4.3) 

where ki^ and /c2/i denote the momenta of the incoming electron and positron, = ki^ + 
k2^, and the contributions 

to the hadron tensor. They will be given below. Performing the tensor contractions in 
d = 4 — 2e dimensions, one obtains 

^«i4;v,=6^e^(l-6)5]H«g. (4.5) 

X 

The contribution to (4.5) from the QQqq final state is closely related to the integrated 
antenna function B^^-q given in eq. (3.19). After restoring all couplings and color factors 
we find 

a^Nf g2(3_2e) (87r2(4^)-^e7B^)2 QggQvy.y.F, v ; 

where Cp = (iVc ~ ^)/{'^^c) and Tr = i. The expressions for \Mqq\'^ and P2{q^,m) are 
given in (2.4) and (3.3), respectively. The second normalization factor in (4.6) is obtained 
in straightforward fashion; it reflects the relation between the decay rate of a virtual photon 
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and the integrated antenna function, see (2.5), (3.1) and (3.2). In our calculation of the 
antenna function the hadronic tensor (4.4) was contracted with — instead of L^y. This 
is corrected by the additional factor l/(q^(3 — 2e)). 

The two-particle contribution n*^*^ can be expressed by the Dirac and Pauli heavy 
quark form form factors Fi and F2- 



a? AT 



/ 3 - 2e 



(1 + y2 + 2/(10 - 8e)) , „ ^ .... 



(4.7) 



The contributions required in (4.7) can be read off from the expressions for the UV- 
renormalized form factors above threshold given in [50]. 

In [50] the renormalization constants for the heavy quark mass and wave function were 
defined in the on-shell scheme, whereas the renormalization of the strong coupling constant 
and the gluon wave-function was performed in the MS scheme. In order to obtain an LSZ 
residue equal to one, we apply on-shell renormalization for the external gluon, too. (Nom- 
inally, this avoids contributions from three-particle cuts.) This change of the renormal- 
ization scheme as compared to [50] leaves the two-loop contributions oc a'^Nf unchanged. 
However, it changes the QQg renormalization constant ZiF{e, fi'^ /m?) as compared to the 
one of [50] at the one-loop level by the additional term 



r (1 + e) . (4.8) 



This change induces a counterterm contribution proportional to a^Nj from the three par- 
ticle QQg final state, which reads: 

where -^QgQ is the integrated massive tree-level three parton quark-antiquark antenna as 



given, e.g., m 

Adding the contributions (4.6), (4.7), and (4.9) to (4.5), all IR poles cancel. This 
provides a strong check for the IR divergent part of the integrated antenna function given 
in (3.19). 

Next we compare our result for Ra^^a ^ with the result of of [76], which was obtained in 
d = 4 using different techniques. Introducing the QCD coupling ag and the appropriate 
color factor in eq. (42) of ref. [76] the relevant part of this equation becomes 



i W In (^) + + u;(3 - u,^) (/f ^ + ) + 



3^(0) 
2 



, (4.10) 



where 

w = l^ = VT^, (4.11) 
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and the explicit expressions of the functions W, fi^\ and f2^^ can be found in reference 
[76] . In [76] the result for f^^^ is expressed in terms of the integrals 



Jo 



^ arctan(^j;) 



3.2 _|_ ^2 



niv, = I dx ^"^^ (4.12) 



niv, e, x) = / dx 

Jo 



J.2 _j_ j^2 

1 ^„l^2 I c2 



ln(x + ^ ) arctan(x2;) 



J.2 _|_ ^2 



In order to be able to compare our result with (4.10) in analytic fashion, we have computed 
the particular combination 
of polylogarithms. We find 



the particular combination of these integrals, which appears in the function in terms 



-r3(l, 0, w) + T;[l, -,w\ -tJi,w, -\ +2\n{w)T2{l,w)-^T*(l, - 
^ w J \ w J 2 \ w 

Lis! ) -Lisf ) +ln(u;) (Lisf ) - Lis 



1 — wJ \l + wj \ \1 + w J \ 1 — w 

1 o 1 o 7r2 7r2 

+ - ln'^(l + w)- - ln%l -w)-— ln(l + w) - — ln(l - w) 
6 6 3 6 

+ — ln{w) + - In^iw) (ln(l - w) - ln(l + w)) + ttG , (4.13) 

where G is Catalan's constant. With this formula we find agreement^ between our result 
(4.5) and the result (4.10) of [76]. 



5 Summary and Outlook 

As a first step towards extending the antenna subtraction method to NNLO QCD reac- 
tions with massive quarks, we have determined the real radiation antenna function and its 
integrated counterpart for reactions of the type S — QQqq, where 5 denotes an uncolored 
initial state. We were able to determine the integrated antenna function in completely an- 
alytic fashion, namely in terms of harmonic polylogarithms, for which efficient evaluation 
codes are available. We checked our results by computing the contribution proportional 
to a'^e'qNf to the inclusive heavy-quark pair production cross section in e'^e~ annihilation 
via a virtual photon and by comparison with results in the literature. 

An obvious next step in this line of investigation is the determination of the antenna 
function and its integrated version for S — t- QQgg. The results of this paper indicate that 
for the QQgg final state, the integrated antenna function can also be obtained analytically 
in a relatively compact form. 

^In the course of this comparison, we found that eq. (23) of [76] contains a typographical error. In the 
fifth line of this equation, Inp^ must be replaced by In^p. 
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A The master integrals 

In this appendix we give analytic results for the five master integrals Ti, T2, T3, T4, T5 to 
the required orders in e. The integrals are given in terms of the variable 



(A.1) 



1 + ^1-^ 

The harmonic polylogarithms are given in the notation of [68]. The e-expansion is needed 
to order for Ti, T2 and Ts and to order e for T4 and T5. For convenience, we define 

T,{q^y,^.^e) = C^e) (^) ^ ^n j(i + ,)6 (A-2) 



with 



We find 



Cie)=[^y. (A.3) 



Ti(y, e) = -^(1 + y) (-1 - 23y + (-34 + Att^) / + (34 + 4tt^) + 23/ + 
+ y (1 + 5y + 6/ + 5^^ + /) H{0; y) - Ay^{l + y f H{-1, 0; y) 
+ 2y^{l + yfH{0,0;y)} 

+ e|4y2(l + yf H{-3; y) + 2y (l + Sy + Qy^ + 5^^ + y^) H{-2; y) 

- i(l + y) (-1 - 23y - 2 (l7 + Svr^) y^ + (34 - 16^^) ^3 ^ 23y4 + y^) H{-1; y) 

- ^y (-45 + (-183 + 8^2) ^ + 2 (-35 + Stt^) y^ + (45 + 8^^) ^3 

+51y4 + 4y5) H{0; y) - ^ (-1 - 24y - 57y2 + 57y^ + 24y^ + y^) H{1; y) 

+ lOy (1 + 5y + 6y2 + Sy^ + y^) H{2- y) + 20y2(l + yf H{3; y) 

+ 28y2(l + yf H{-2, 0; y) - 8y\l + yf H{-1, -2; y) 

+ 2y (5 + 17y + 16y2 + 17y3 + 5y^) H{-1, 0; y) - 40y2(l + yf H{-1, 2; y) 

- y(l + yf (1 - 5y + y^) H{0, 0; y) - 56y2(l + yf H{-1, -1, 0; y) 
+ 12y2(l + yf H{-1, 0, 0; y) - 6y2(l + yf H{0, 0, 0; y) 
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+ ^ (71 - 24 (-71 + Stt^) y - 24 (71 + Stt^) / - 71^^ 
+ (4047 - 4567r^ - 720 C(3)) - 120y^ (Stt^ + 12C(3)) 
-3/ (1349 + 1527r2 + 240 C(3)) ) | 

+ - ^ (-445 + Qtt^) + A (_623 + Ihir'') + (-12^^ _ 2Ay^ 

-12y4) i?(-4; y) + {-2y + Gy^ + 16/ + - 2y5) F(_3; y) 

+ (^15y - ^ (-183 + SOtt^) / _ ^ (_7 + i6;r2) - ^ (9 + IGtt^) 

-17y5-^) i^(-2;y) 

/ 355 710y 6745y^ 6745y-^ 710^^ 355y6 \ 
V'se" 3 12 3 36~ ) ^ ' 

+ (75y + 305y^ + ^ _ 75y^ - SSy^ - ^) H{2- y) 

+ (-lOy + 30y2 + SOy^ + 30y^ - lOy^) H{?>- y) 

+ (-60y2 - 120y^ - 60y^) H{^-y) 

+ (8y2 + 16y3 + 8y4) iJ(-3,-l;y) 

+ (-132y2 - 264y3 - 132y^) iJ(-3,0;y) 

+ (40y2 + 80y3 + 40y^) H{-Z, 1; y) 

+ (56y2 + 112y3 + 56y^) H{-2, -2; y) 

+ (4y + 20y2 + 24y3 + 20y^ + Ay^) H{-2, -1; y) 

+ (-38y - 78y2 - 32y3 - 78y^ - 38y5) H{-2, 0; y) 

+ (20y + 100y2 + 120y3 + IQOy^ + 20y^) H{-2, 1; y) 

+ (280y2 + 560y3 + 280y^) H{-2, 2; y) 

+ (24y2 + 48y3 + 24y^) H{-1, -3; y) 

+ (20y + 68y2 + 64y3 + 68y^ + 20y^) H{-1, -2; y) 

/I 1 , o 3207r2y3 1 , , 
+ f 3 + 8y + - (57 + leOvr^) y^ + ^ + - (-57 + IGOtt^) y^ 

-Sy'-^) H{-l,-l-y) 

+ (^-^ + 19y - ^ (-408 + 77r2) y^ - ^ (-136 + y^ - ^ (-180 + Itt'') y^ 
-13y5-^) i?(-l,0;y) 

+ Q + 40y + 95y2 - 95y4 - 40y5 - H{-1, 1; y) 

+ (lOOy + 340y2 + 320y3 + 340y^ + lOOy^) H{-1, 2; y) 
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+ (1202/2 + 240y3 + uOy^) H{-1, 3; y) 

+ ^(0,0; y) 

+ Q + 40?/ + 952/2 _ 95/ - 402/5 - ^) (1, -1; y) 

/lO 10u^\ 
+ ^ + 80y + 190y2 - IQO/ _ go/ _ ^ O; y) 



3 " " " " 3 , 

+ + 200y + 475/ - 475/ - 200/ - H{1, 1; y) 

+ {20y + 100/ + 120/ + 100/ + 20/) H{2, -1; y) 
+ (40y + 200y2 + 240/ + 200/ + 40^^) H{2, 0; y) 
+ (lOOy + 500y2 + 600y^ + 500/ + lOOy^) H{2, 1; y) 
+ (40y2 + 80y^ + 40/) H{3, -1; y) 

+ (80y2 + 160/ + 80/) H{3, 0; y) + (200^^ + 400/ + 200/) H{3, 1; y) 

+ (296/ + 592/ + 296/) H{-2, -1, 0; y) 

+ (-36/ - 72/ - 36/) H{-2, 0, 0; y) 

+ (-16/ - 32/ - 16/) H{-1, -2, -1; / 

+ (264/ + 528/ + 264/) H{-1, -2, 0; y) 

+ (-80/ - 160/ - 80/) H{-1, -2, 1; y) 

+ (-112/ - 224/ - 112/) H{-1, -1, -2; / 

+ (92y + 236/ + 160/ + 236/ + 92/) H{-1, -1, 0; y) 

+ (-560/ - 1120/ - 560/) H{-1, -1, 2; y) 

+ {-6y + 18/ + 48/ + 18/ - 6/) H{-1, 0, 0; y) 

+ (-80/ - 160/ - 80/) H{-1, 2, -1; y) 

+ (-160/ - 320/ - 160/) H{-1, 2, 0; y) 

+ (-400/ - 800/ - 400/) H{-1, 2, 1; y) 

+ {y- 19/ - 36/ - 19/ + /) H{0, 0, 0; y) 

+ (-592/ - 1184/ - 592/) H{-1, -1, -1, 0; / 

+ (72/ + 144/ + 72/) H{-1, -1, 0, 0; y) 

+ (-28/ - 56/ - 28/) H{-1, 0, 0, 0; y) 

+ (14/ + 28/ + 14/) H{0, 0, 0, 0; y) 

(1\ 2 , 9, 2 , a, . 71/ 

+ (^3g - 3 (-71 + levr^) y - - (71 + levr^) / - ^ 



36 

- ^y' (-1349 + 5127r2 - 384 C(3)) 
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-^2/4 (1349 + 5127r2 - 384C(3)) - ^y^ (ITtt^ - 24C(3))) Hi-l;y) 

(3 1 llv^ 1 

-- (-47 + 7r^)y-- (1849 + gTr^) - - (_n2i + 10977^ + 192 C(3)) 

-^y^ (4275 + 1097r2 + 192 C(3)) - (3393 + 2497r2 + 576 C(3))j H{0;y) 



+ ^59185 _ 602l7r2 + TGtt^ - 18720 C(3)) 

+ (-59185 + 230l7r2 + 767r^ - 18720 C(3)) 

+ ^y^ (-2497r2 + 197r^ - 3204 C(3)) 

lo 

+ (_3i;l5 ^ 2407r2 - 324 C(3)) 

- ^y (-3115 + 1987r2 + 324 C(3)) |. (A.4) 

Uq^y, /X, e) = C\e) ,nj(i^^^)8 T-^Cy. ^) > (A.5) 

with 

r2(y, e) = -^(1 + y) (-1 - 232/ + 67y2 + (149 - 24^^^ _ (^^49 ^ 247r2) / - 67^^ 
+23y6 + y7) + 1 _ i3y3 _ 14^4 _ -^3^5 ^ ^7) ^^q. 
+ 4y3(l + y)2 0; y) - 2y\l + yf H{0, 0; y) 

+ e| - 4y3(l + y)2 //(_3; y) + 1 (y - 132/3 _ 142/^ _ 132/^ + y^) Hi-2; y) 

- ^(1 + 2/) (-1 - 23y + 67y^ + (149 + 967r2) 2/^ + (-149 + 967r2) - 67y^ 
+23y6 + /) i?(-l;y) 

+ ^y (45 - 30y + (-545 + 487r2) y^ + (-6 + 96tt^) y^ + (319 + 487r2) y^ 
+U6y^ - 51y^ - 4y^) H{0;y) 

- ^ (-1 - 24y + 44y2 + 216y3 - 216y5 - 44y6 + 24y^ + y^) H{1; y) 

+ \{y- ^^y^ - 1 V - 13y' + y^) H{2- y) - 20y3(l + yf H{3; y) 

- 2^y\l + yf H{-2, 0; y) + 8y3(l + yf H{~1, -2; y) 

+ ^y (5 - 23y2 + 4y3 - 23y^ + Sy^) H{-1, 0; y) + 40y3(l + yf H{-1, 2; y) 

- ^2/(1 + y)' (1 - 2y + 32y2 - 2y3 + y^) i/(0, 0; y) + 56y=^(l + yf H{-1, -1, 0; y) 



12y3(l + y)2 i7(-l, 0, 0; y) + 6y=^(l + yf H{0, 0, 0; y) 
+ ^ (71 - 24 (-71 + 3tt^) y - 2704y^ + 2704y^ - 24 (71 + Stt'^) y^ - 71y^ 
+ 24y3 (-601 + eOvr^ + 180 C(3)) + 24y5 (601 + 607r^ + 180 C(3)) 
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+24/ (797r2 + 360C(3))) | 

+ I I2y\l + yf H{-4; y) - ^y{l + yf {l - 2y + 32/ _ 2/ + /) H{-2>- y) 

+ ^y (45 - 30y + 5 (-109 + 96tt^) / + (-6 + 9607r^) / + (319 + 4807r2) / 

18 

+146/ - 51/ - 4/) i?(-2;/ 

- ^ (-71 - 1704y + 2704y^ + 14424y^ - 14424^^^ - 2704y*^ + 1704y'^ + 71y^) H{1; y) 

- ^y (-45 + 30y + 545^^ + Qy^ - 319/ - 146/ + 51/ + 4y^) H{2; y) 

lo 

-\{y + 29/ + 60/ + 292/5 + ^(3. y) + 60/(1 + yf H{A- y) 

- 8/(1 + yf H{-3, -1; y) + 132/(1 + yf H{-3, 0; y) 

- 40/(1 + yf H{-3, l;y) - 56/(1 + yf H{-2, -2; / 

+ l{y- 13/ - 14/ - 13/ + /) H{-2, -1; y) 

- (19 + 47y2 + 252/ + 47/ + 19/) H{-2, 0; y) 

+ y (y - 13/ - 14/ - 13/ + y7) H{-2, 1; y) 

- 280/(1 + yf H{-2, 2; y) - 24y3(l + yf H{-1, -3; y) 
+ ^y (5 - 23/ + 4y3 - 23y^ + 5y6) H{-1, -2; y) 

- -^(1 + y) (-1 - 23y + 67y2 + (l49 + 96071^) y^ + (-149 + 96071^) y^ 
18 

-67/ + 23/ + y^) F(-l,-l;y) 

+ (-7 + 57y + 182y2 + 7 (-85 + Gtt^) y^ + 6 (87 + Mtt^) y^ 
lo 

+ (269 + 427r2) y^ + 358y6 - 39y^ - lly^) H{-1, 0; y) 

- ^ (-1 - 24y + 44y2 + 216y3 - 216y5 - 44y6 + 24y^ + y^) H{-1, 1; y) 

+ y y (5 - 23y2 + 4y3 - 23y^ + 5y6) H{-1, 2; y) - 120y3(l + yf H{-1, 3; y) 

+ ^y(l + / (-45 + 51y - 2 (62 + 2l7r2) / - 2 (211 + 2l7r2) / 
-83/ + / + 2y6) H{0,0;y) 

- 4 (-1 - 24y + 44y2 + 216y3 - 216y5 - 44y6 + 24y^ + y^) H{1, -1; y) 
lo 

- ^ (-1 - 24y + 44y2 + 216/ - 216/ - 44y6 + 24y^ + /) H{1, 0; y) 
25 

- — (-1 - 24y + 44y2 + 216y3 _ 216/ - 44y6 + 24y^ + y^) H{1, 1; y) 
+ f (2/ - 13/ - 14/ - 13/ + y') H{2, -1; y) 
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+ Y (y - 13?/' - 14/ - 13y5 + y') H{2, 0; y) 

+ - (y - 13y3 _ 1 V - 13?/5 + /) H{2, 1; y) - 40y3(l + yf i/(3, -1; y) 

- 80y3(l + y)2 /7(3, 0; y) - 200^^(1 + y f H{3, 1; y) 

- 2962/=^(l + i7(-2, -1, 0; y) + 36y\l + yf H{-2, 0, 0; y) 

+ 16y3(l + y)2 _2, -1; y) - 264y3(l + y)^ H{-1, -2, 0; y) 

+ 80y3(l + y)2 -2, 1; y) + 112y3(l + y)^ H{-1, -1, -2; y) 

+ ^y (23 - 5y2 + 196y3 - + 23y6) -1, 0; y) + 560y=^(l + yf H{-1, -1, 2; y) 

- y(l + yf (1 - 2y + 32y2 _ 2y3 + q, 0; y) + 80y^{l + yf H{-1, 2, -1; y) 
+ my^{l + yfH{-l,2,0;y) 

+ 400y=^(l + yf H{-1, 2, 1; y) + 1 (y + 113y3 + 208/ + 113/ + /) H{0, 0, 0; y) 

+ 592y=^(l + yf H{-1, -1,-1, 0; y) - 72y=^(l + yf H{-1, -1, 0, 0; y) 
+ 28y=^(l + yf H{-1, 0, 0, 0; y) - 14y=^(l + yf H{0, 0, 0, 0; y) 

+ ^ (71 - 24 (-71 + 167r^) y - 2704y2 + 2704y^ - 24 (71 + levr^) y^ - 71y^ 

+ 24?/ (-601 + 12471^ - 288 C(3)) + 24y^ (601 + 1247r2 - 288 C(3)) 

+96/ (1977^ - 144 C(3))) H{-l;y) + ^y (l269 - 2550y + 8266/ 

- 5547y^ - 284y^ - Stt^ (9 - 453y2 - 718y^ - 453y^ + 9y^) 
+ 6y3(3535 + 1152 C(3)) + y2(-15937 + 3456 C(3)) 
+/(41759 + 3456C(3))) i?(0;y) 

+ -i- (3115 - 105152/ + 105152y6 - 3115/ - 13687r^y^(l + yf 

+ 37r^ (-21 - 1584y + 1740y2 + 20088y^ + 2352y^ - 11448y^ - 4684y6 
+1920/ + 125y^) + y(74760 - 7776 C(3)) + 268704y^ C(3) 

- 24y^(3115 + 324 C(3)) + 24y3(-25145 + 7992 C(3)) 

+24y^(25145 + 7992C(3))) |. (A.6) 
Tsie, y, e) = C\e) ,njfi + ^)8 r.{y, e) , (A.7) 

with 

TsCy, e) = -^(1 + y) (-3 - 109y + (-443 + 247r2) / + (-277 + 727r2) y^ 
+ (277 + 727r2) y^ + (443 + 247r2) y^ + lODy^ + 3y^) 
+ ^y (2 + 18y + 42y2 + 57/ + 42/ + ISy^ + 2y6) H{{)- y) 

- 4/(1 + yf H{-1, 0; y) + 2/(1 + yf H{0, 0; y) 
+ e\4y\l + yfHi-3;y) 
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+ -y (2 + 18y + A2y^ + 57y3 + A2y^ + 18/ + 2y^) H{-2; y) 

+ ^ (3 + 112y + 24 (23 + A-k"^) y^ + 48 (15 + Stt^) y^ + STBTr^y^ 
+48 (-15 + 8tt^) + 24 (-23 + iir^) / - 112/ - 3/) H{-1; y) 
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